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SUMMARY 

A creep damage accumulation model is presented that makes use of the 
Kachanov damage rate concept with a provision accounting for damage that 
results from a variable stress history. This is accomplished through the 
introduction of an additional term in the Kachanov rate equation that is 
linear in the stress rate. Specification of the material functions and param- 
eters in the model requires two types of tests constituting a data base: 

(1) standard constant-stress creep rupture tests, and 

(2) a sequence of two-step stress creep rupture tests 

INTRODUCTION 

One of the primary failure modes considered in the ASME Code Case N-47 
(Ref. 1) for structural components in elevated-temperature service is that of 
creep rupture under quasi-steady, long-term loading. The Code Case specifies 
that predictions of time to failure under such loading conditions should be 
based on the "time-fraction law," originally proposed by Robinson (Ref. 2), 



in which t^ denotes the time to failure and t R (o) represents the "creep- 
rupture curve" determined from uniaxial tensile creep tests at constant stress. 

Several authors (e.g.. Refs. 3 to 5) have pointed out that Eq. (1) is not 
strictly satisfied for many structural alloys under variable stress conditions. 
For example, it is oftet observed that 


( 2 ) 



in uniaxial two-step creep rupture tests involving a step-up in tensile stress 
and 



for a stress down-step*. 

Here, a creep damage model is proposed similar to that originally pro- 
posed by Kachanov (Ref. 6) - and employed or extended by numerous workers 
(e.g.. Refs. 7 to 9) - but which includes an additional term depending linearly 
on the stress rate. This resjlts in a time-independent contribution to the 
predicted damage arising from variations in stress in an already creep damaged 
material. The model is stated here for isothermal and noncyclic conditions; 
extensions to nonisothermal conditions and to reverse stressing will be sub- 
jects of continued research. 

The proposed damage rate equation is shown to be equivalent to a damage 
law of the form 



in which t^ and t^(o) have the same meanings as in Eq. (1) and o is a 
functional of the stress history o(t). For some stress histories o is nega- 
tive and Eq. (4) predicts a comparatively shorter time to failure than Eq. (1); 
in other cases a is positive and Eq. (4) predicts a longer time to failure. 


♦Notably, this effect of stress sequencing is opposite to that ordinarily ob- 
served in fatigue testing. 
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The implementation of the model requires a data base made up of two types 
of tests: 

(1) standard constant-stress creep rupture tests, and 

(2) variable (two-step) stress creep rupture tests. 

A unique feature here is that the variable stress experiments are used 
not solely as tests to furnish evidence for verification of the model, as is 
commonly the case, but instead they make up an essential part of the data base 
for establishing the functional forms and parameters of the model. 

As a background for the development of the proposed model, the connection 
between the Kachanov damage model and the time-fraction law, Eq. (1), will 
first be demonstrated. The proposed model will then be stated along with a 
discussion of the testing required for its specification. 

THE KACHANOV MODEL AND THE TIME-FRACTION LAW 


The Kachanov damage model, giving the rate of degradation of a material 
under creep at elevated temperature, is frequently stated in the separable form 


f(q) 

g'U) 


(5) 


in which is the material continuity and o is the applied stress. For a 
material completely intact <j» = 1, and for complete destruction ^ = 0. 

The rupture life t^ under constant stress is found by integration of 
Eq. (5) and is 


t R (o) 


g(l) - g(o) const. 
?(o) = T^T 


( 6 ) 


In the special case where f = Co n , Eq. (6) takes the familiar form 

^ , . const. 

M 0 ) = n 


The failure time t^ under an arbitrarily varying stress is found from 
Eqs. (5) and (6) as 


(7) 
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which is identical to the time-fraction law, Eq. (1). Whatever the functions 
f and g, Eq. (5) demands the satisfaction of Eq. (1) at failure. The 
Kachanov creep damage model, expressed in the separable form of Eq. (5), 
cannot lead to the conditions stated in Eqs. (2) or (3). 

THE PROPOSED DAMAGE MODEL 

Here, a damage rate equation similar in spirit to that of Eq. (5) is 

introduced that contains an additional term proportional to the stress rate 
That is, we propose the rate equation 

, t n 

= h(o,i|))o - C (10) 

in which the second term is taken to be a special form of the Kachanov equa- 
tion, Eq. (5), whereas the first term is taken linear in the stress rate a. 
The first term represents a time-i ndependent contribution which arises from 
changes in applied stress; the second term represents the usual time-dependent 
contribution as proposed by Kachanov. The constants C, n, and m and the 
function h are determined by experiment as described below. 

Equation (10) is consistent with the view that creep damage occurs as the 
result of the creation of interior voids or cavities. In the context of step- 
wise creep tests, the first term in Eq. (10) allows for either further damage 
or healing of prior damage as a consequence of changes in stress. This can be 
thought of as resulting from further wid ening or partial closing of existing 
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voids depending on the direction of the stress change. So that stress changes 
do not effect the yet undamaged material, it is reasonable to restrict the 
function h such that 


h(o,l) = 0 (11) 

With this restriction, it is seen that in ordinary creep rupture tests, 
in which an undamaged specimen is brought relatively abruptly to a constant 

stress and held (i.e., a = 0 ), there is no contribution from the first term in 

Eq. (10) throughout the test and we can write 

» 

0 

# K 

( 12 ) 
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Co n dt 


or 


tn ( o ) - 


(13) 


C(m ♦ l)o 

Optimal values of l/[C(m ♦ 1)] and n can thus be determined using Eq. (13) 
and data pairs in the form (t^, o) from standard creep rupture tests. 

The variation of the material continuity 4 ) in creep rupture tests is 
likewise found by integrating Eq. (10) with the first term absent. Thus 

\1/ (m+1 ) 




(14) 


Only for the very special case m = 0 , corresponding to the second term 
in Eq. (10) being independent of 41 , is the diminuation of material continuity 
(or the accumulation of damage D = 1 — u» ) linear in time. With m > 0 , 41 
varies nonlinearly in time, slowly at first but at a higher rate as internal 
damage occurs. 

Hypothetical trajectories of two creep rupture tests at stress levels of 

plot of Fig. (1) . 


oj and o,, are shewn as the solid lines in the 




Each curve originates at the intersection of = 1 and / = o and each 

terminates, at failure, in the plane = o. The locus of rupture points in 


<)i = o is defined by 


/ 


dt 


= 1 . 


We now assume that the stress dependence of the minimum creep rate cm 
in constant-stress creep tests can be adequately represented by the Norton law 


(Ref. 10) 


£ = Ao 
m 


(15) 


where A and N are known constants. The subsequent development is not de- 
pendent on the particular choice of the Norton law; it has been chosen because 
it has wide application and is simple. We further assume that the subsequent 
acceleration of creep in such tests is due to internal damage accumulation* 
and that, following Kachanov, Eq. (15) remains valid with o replaced by the 
"effective" stress o/ 4 >, i.e., we take 


/ \N 


(16) 


Dividing Eq. (15) by Eq. (16) and solving for 4 * , we get 

1/N 



(17) 


indicating that we can assign values of 41 along a given creep-rupture curve. 


The minimum creep rate c m can be readily identified, and c can be deter- 
mined at each instant of time along the curve. 

Eliminating 41 from Eqs. (14) and (17) we get: 


*If, as is usually the case, the creep rupture tests are conducted under 
constant load rather than constant stress, the contribution of the creep 
acceleration related to geometry change must be accounted for. 
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m 




N/(m+l) 


( 18 ) 


As (e m /e) is measurable along a given creep curve, we can generate data sets in 
the form (e m /e, . , t) directly from creep-rupture tests and use these with 
Eq. (18) to obtain optimal values of N/(m + 1 ). with N known, this estab- 
lishes m; with m and l/[C(m + 1 )] in Eq. (13) known, C is determined. 

Thus, n, m, and C are known and the second term in Eq. (10) is fully speci- 

fied from creep-rupture data alone. 

Before proceeding with the determination of h ( 0 , 4 .) in Eq. (10), we 
consider the predicted response of a two-step stress test when h = 0 and 
Eq. (10) is of the classical Kachanov form. In this case, the image of a 
typical two-step creep rupture test is shown as OABC in Fig. 1 . As the abrupt 
step from oj to o 2 occurs (AB), 4 . remains constant. Note that the pro- 
jections of the step test OABC as well as those of the constant stress tests at 

oj and are identical in the 
constant or variable stress, 


•■/v 


plane. Moreover, for each test. 



(19) 


at failure. 

Next, we consider stepped creep rupture tests designed to establish the 
form of the function ( 1 ( 0 , 4 ). The image, predicted by Eq. (10), of a test 
involving a step increase from oj to o ? is indicated as OABC in Fig. 2 . 
The stress is held constant at o^ in 0 < t < t^ (t^ being the time corre- 

sponding to point A). At t = t^ the stress is abruptly increased to o 2 
(along AB) and again held constant. We suppose that creep rupture occurs at 
the time t = t c (corresponding to the point C in Fig. 2). Here, * is re- 


7 


presented as 
abruptly to 


decreasing (indicating further damage) as the stress is increased 
a?. Correspondingly, 



( 20 ) 


at the failure time t^ = tc. 

The case where a step down in stress is made from a 2 to o 1 at t = t A 
is depicted in Fig. 3. Here, * is represented as increasing (indicating 
partial healing) as the stress is decreased and, 



(21) 


at the failure time t f = tc. 

Over the intervals when the state point follows the path segment OA or BC 
in Figs. 2 or 3, o = 0 and the first term in Eq. (10) does not contribute to 

During initial load-up * = 1, corresponding to the undamaged material 
condition, and likewise the first term does not contribute, in accordance with 
Eq. (11). However, as the abrupt stress change AB is made, the material has 
incurred prior damage and, as the time interval of loading is short, the first 
term in Eq. (10) dominates; the continuity *, then increases or decreases, 
depending on the nature of h(a,ii>) and the stress history. 

In light of the foregoing discussion, * A , (e.g., in Fig. 2), is found by 
integrating Eq. (10), including only the second term, thus 



or 
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( 23 ) 


/ t. \ 1/(m,1) 

*A “ (l - Wj 

Sim.iarly, 4 «b is found by 



The change in 41 during the abrupt step change from to o 2 is thus the 
difference of Eq. (23) and Eq. (25). Equivalently, we know the points A' 

(oj,t|/ A ) and B' (o 2 , 4 > B ) as projected into the 0,41 plane (Fig. 4). 

If we were to perform a sequence of step tests of this kind, each starting 
at o = oj and with a step (up or down) at t = t^, we would map out a curve 
in the a,* plane containing A' and B'. Further, if we conducted several 
sequences of step tests commencing at other stress levels and with stress 
changes at other times, we would eventually map out an entire family of curves 
in a region of the 0 , 4 / plane (Fig. 5). Figure 5 is schematic and real data 
would likely show considerable scatter, nevertheless, it is reasonable to ex- 
pect that the underlying trends can be satisfactorily represented by a family 
of curves 

p(o, 4 > ) ” const. (26) 

which fits the data in a least-squares sense. 

As the abrupt stress change occurs in any of the step-stress tests, the 
state point moves along one of the curves represented by Eq. (26). As discussed 
earlier, the first term in the rate Eq. (10) governs, i.e., we have 

41 = h(o, 4 >)a (27) 
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We can eliminate time in Eq. (27) and write 

$7 n (28) 
Now, we want the experimental curves represented in Eq. (26) to be the 
integral-curves of Eq. (28), that is to say, we want the function hfg,*) to 
be given by 



h(o,n») 


(29) 


With p(o, 4 «) specified, and thus ap/aa and ap/a* known, we then have 
an explicit form for h( a ,^ ) . 

This completes the specification of the rate equation, Eq. (10), and we 
can supposedly use it to predict the time to creep failure under an arbitrary 
stress history a(t). In general, this is done by numerically integrating 
Eq. (10) to determine the time at which * » o. 

The approach followed here in determining the function h ( a , 4 » ) on the 
basis of two-step creep rupture tests is analogous to that followed by Leckie 
and Ponter (Ref. 11) in obtaining representations of creep deformation under 
variable stress conditions. In their case, stepped creep tests provided a 
description of how to move from one constant stress creep curve to another as 
the stress varied. Here, stepped creep rupture tests similarly provide a 
description of how to move from one constant stress rupture curve, in the 



space of Fig. 2 or 3, 


to another as the stress is varied. 


The success of the model proposed here depends on the ability to conduct 
a sequence of accurate and repeatable two-step stress rupture tests. Assuming 


this can be done, the model must, at the very least, give back accurate pre- 
dictions of creep rupture under stepped loading, since such tests comprise its 
data base. In principle, the model should also provide good failure predic- 
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tions under general (tensile) variable stress conditions, although the proof 
of this must obviously come from experiment. 

THE PROPOSED MODEL AND THE TIME-FRACTION LAW 
To demonstrate that Eq. (10) is equivalent to the form shown in Eq. (4), 
we integrate Eq. (10) for an arbitrary stress history o(t) as follows: 



in which 

r°< 

a = (in + 1) I 4 i IT \i( o , 4 » )do (33) 

A 

Note that the integral in Eq. (33) is path dependent so that the value of o 
depends on the stress history o(t). 

SOME HYPOTHETICAL FORMS OF p(o,*) 

The key to the behavior under variable stress lies in the form of p(o,*) 
in Eq. (26). Experiments are likely to show that the p curves of Fig. 5 are 
characteristical ly different in different regions of the o,^ plane, reflect- 
ing distinct rupture mechanisms. In any case, if a sufficient number of vari- 
able stress rupture tests are conducted spanning the relevant regions of the 
a,<|i plane, this information should be inherently built-in to the model. 


In the absence ot sufficient variable stress tests, it is instructive to 
hypothesize some simple forms of p(o,i|») and examine the consequences on the 
predicted creep rupture behavior. 

In regions of the a,i^ plane where Eq. (26) represents a family of 
horizontal lines (Pig. 6), h = o in Eq. (10) and the model is equivalent to 
the usual t ime-fraction law expressed in Eq. (1). 

An interesting case is provided by the family of curves 


P(o.l>) = — = z~ = 


const. 


in which s and a Q are constants. We then have 


ap _ 1 - ii»_' 

ao 2 


i£ = - i .s-l 


so that 


Mo, * ) ■•kfet-) 


Note that, in agreement with Eq. (11), 


h (o,l ) = o 


Equation (10) thus becomes 


• i /i - „ s y r „ n 

* = ~ 17 l"rr ) ° • c ns 
\ * / * 


The family of curves corresponding to s = 2 and several values of o ir 

o 

Eq. (34) is shown in Fig. 7. Using these and appropriate values for m, a., 
and o f , we can calculate the time-fraction integral, Eq. (32), for any given 
stress history. Consider, for example, the two-step tests depicted in Fig. 8; 
the first involves a step-up in stress from a- = 140 to a f = 175 MPa (Fig. 



8(uji) and the second a step-down from = 140 to = 105 MPa (Fig. 8(b)). 
With m taken as m = 1.5, values of the time-fraction integral have been cal- 
culated and are given in Tables I and II for the stress histories of Figs. 8(a) 
and (b), respectively. The quantity a in the tables denotes the fraction of 
the rupture life t^(o.) = t R (140) at which the step up (or down) in stress 
is made. When the step is made relatively early in life (e.g., a = 0.2), the 
value of the integral at failure is not very different from unity; when the 
step is made relatively late in life (e.g., e = 0.8), the integral may differ 
substantially from unity. 

The point 0 in Fig. 7 represents the point in the a, 4 » plane from which 
the abrupt stress step is made corresponding to b = 0.2 in Table I or II. As 
the stress is increased from 140 to 175 MPa, the state point moves from 0 to a. 
As the stress is decreased to 105 MPa, the point moves to b. Only slight addi- 
tional damage or healing is incurred during the stress change. The predicted 
values for the time-fraction integral at failure are 



for the step-up in stress, and 



for the down-step. 

The point O' in Fig. 7 indicates the point from which a step is made with 
b - 0.8. The material has now incurred considerable prior creep damage, and a 
step-up in stress to 175 MPa moves the state point to a', indicating substan- 
tial additional damage as the stress is increased. The value of the time- 
fraction at failure in this case is 
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f: 


dt 


t7(oT 


= 0.85 


(42) 


Note that here, following the expenditure of 0.8 of the rupture life at 
o = 140 MPa, a stress increase to approximately 193 MPa or greater causes 
failure immediately as the stress is changed, corresponding to * » 0 . 

A step down in stress from O' moves the state point to b' providing par- 
tial healing, and the integral at failure is 

’t. 


/ 


dt 

tTftT 


= 1.18 


(43) 


Additional forms of p(o,*) can be hypothesized and calculations made 
s^ilar to those above, however, such exercises would probably not be fruitful 
at present. It is preferable to wait until sufficient experimental data can 
be generated and the actual forms of p(o,„,) and h(o,*) determined. 

It is possible that a single family of curves in the o,* plane is not 
sufficient to represent the effect of both stress increases and stress de- 
creases on the creep damaged material. If this were the case, the same prin- 
ciples and procedures outlined above are still applicable with the following 

modifications and reinterpretations. The damage rate equation, Eq. (10), could 
then be modified to 


• m 

* = h 1 (o,i|»)<o> + h 2 (o,<j.)<-o> - 


n 

c£ sr 


(44) 


with 

and 

in which 


<x> = x; x > o 
<x> s o; x < o 
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h 1 (o,Y) = - 


( 45 ) 



The family of curves represented by ( 0 *4* ) = const, in the 0,41 plane would 

be determined from a sequence of stepwise stress rupture tests involving a 
step-up in stress and the curves P2(o,<i<) = const, from a sequence with a 
step-down in stress. 

CONCLUSIONS 

The proposed damage model is limited in applicability to a narrow range 
of conditions, i.e., uniaxial, isothermal, variable tensile stress histories. 
Nevertheless, the existing technology supporting the design of elevated tem- 
perature components (e.g., the guidance provided by ASME code case N- 47 ) does 
not provide an adequate methodology for accurately predicting creep rupture 
even under these restricted conditions. 

The attempt here is to provide a simple model which has the potential for 
predicting failure in, perhaps, the most fundamental of problems of high- 
temperature design, creep rupture under quasi-steady, long-term loading. Once 
these predictions can be made with reasonable consistency and accuracy, only 
then does it seem appropriate to concentrate on further complexities such as 
cyclic stressing, multiaxial ity, creep-fatigue interactions, etc. 

As the present creep damage model is based on both standard creep rupture 
tests and stepwise stress rupture tests, it is expected that it should provide 
accurate predictions under general (tensile) variable stress conditions. This, 
of course, must ultimately be demonstrated by experiment. 
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f dt 
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Figure 3 , 


- Representation of step-down creep rupture test 


in 




Figure 4. - Projection of path segment AB 


on a , 4 ) plane 
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